Nonlinear Evolution of the r-Modes in Neutron Stars 
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The evolution of a neutron-star r-mode driven unstable by gravitational radiation (GR) is studied 
here using numerical solutions of the full nonlinear fluid equations. The dimensionless amplitude 
of the mode grows to order unity before strong shocks develop which quickly damp the mode. In 
this simulation the star loses about 40% of its initial angular momentum and 50% of its rotational 
kinetic energy before the mode is damped. The nonlinear evolution causes the fluid to develop 
strong differential rotation which is concentrated near the surface and poles of the star. 
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PACS Numbers: 04.40. Dg, 97.60. Jd, 04.30. Db 

The r-modes of all rotating stars are driven towards 
instability by gravitational radiation (GR) reaction , 
and the strength of this destabilizing force is sufficient 
to dominate over internal dissipation in hot, rapidly ro- 
tating neutron stars || . The growth timescale of the in- 
stability is about 40 s for neutron stars with millisecond 
rotation periods. Thus it is generally expected that GR 
will cause the dimensionless amplitude of the most unsta- 
ble (m = 2) r-mode to grow to order unity within about 
ten minutes of the birth of such a star. The emission 
of GR by this process removes angular momentum and 
rotational kinetic energy from the star. The strength of 
the GR emitted and the timescale for spinning down the 
young neutron star depend critically on the amplitude to 
which the r-mode grows. Initial estimates assumed that 
the amplitude would grow to order unity before some un- 
known process would saturate the mode. With this satu- 
ration amplitude, a neutron star spins down to about one 
tenth its maximum angular velocity in about one year, 
and the GR from this event might be detectable by LIGO 
II W . Stergioulas and Font M find no saturation of the 
r-modes even at large amplitudes in their nonlinear nu- 
merical study. Thus at present no one really knows what 
process will saturate the r-modes, or how large their am- 
plitudes will grow. Here we investigate the growth of the 
r-modes by solving numerically the nonlinear hydrody- 
namic equations driven by GR reaction. 

Neutron stars are compact objects with reasonably 
strong gravitational fields, GM/R ss 0.2c 2 . While these 
objects may contain fluid moving at fairly large velocities, 
v 2 < GM/R, Newtonian theory is expected to describe 
them up to errors of order GM/Rc 2 . For simplicity then 
we study here the nonlinear evolution of the r-modes us- 
ing the Newtonian equations: 

(1) 



dtp + V • ( P v) = 0, 
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where v is the fluid velocity, p and p are the density and 
pressure, <E> is the Newtonian gravitational potential, and 



Fgr is the GR reaction force. The gravitational potential 
is determined by Poisson's equation, 



V 2 $ = 47rGp. 



(3) 



The GR reaction force Fqr due to a time- varying current 
quadrupole (the dominant multipole for the r-modes) can 
be written (see Blanchet || and Rezzolla, et al. 0): 
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where J 22 represents the n'th time derivative of J22, 
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and Y 2 B 2 = 



f x rVl22/v / 6 is the magnetic-type vector 
spherical harmonic. In slowly rotating stars the m = 2 
r-mode projects onto J22 and no other Ji m . The pa- 
rameter k that appears in Eqs. (^) and (|5|) sets the 
strength of the GR reaction force, and has the value 
k = 32y / 7rG/(45v / 5c 7 ) in general relativity theory. For 
practical reasons (see further discussion below), we take 
k to be about 4500 times this value. 

We solve Eqs. (0)~(|^) numerically in a rotating refer- 
ence frame using the computational algorithm developed 
at LSU to study a variety of astrophysical hydrodynamic 
problems Briefly, the code performs an explicit time 
integration of the equations using a finite-difference tech- 
nique that is accurate to second order both in space and 
time, and uses techniques very similar to those of the 
familiar ZEUS code §. 

We find that as written, Eqs. (f|) and (||) for Fgr are 
nearly useless in a numerical evolution. The problem is 
the large number of time derivatives of J22 that appear 
there. In the case of a pure mode with frequency ui this 
problem is easily solved: — (iu)) n J%2- Even when the 
amplitude of the r-mode becomes large, we expect the 
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fluid motion to be dominated by periodic motions at the 
fundamental frequency of the r-mode. Thus, we expect 
the normal-mode expressions for the time derivatives of 
the multipole moments to be reasonably accurate even 
in the nonlinear regime. It is easy to evaluate J22 from 
Eq. @, and can also be expressed as an integral over 
the fluid variables using Eqs. (0) and (0): 



J (1) - 
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Thus we evaluate the time derivatives needed in Eqs. (Q) 
and (||) using = U) A J^ and j!$ = — u e Ji2- We have 
verified that using these approximations our numerical 
code accurately reproduces the analytical description of 
the r-mode evolution and growth due to GR reaction in 
slowly rotating models. 

In order to monitor the nonlinear evolution of an r- 
mode, it will be helpful to introduce nonlinear gener- 
alizations of the amplitude and frequency of the mode. 
Since the r-mode projects primarily onto the multipole 
moment J22, we define the nonlinear amplitude to be 



8^B |J 2 2| 
Q f pr^dV 



(8) 



where VLq is the initial angular velocity and Rq is the 
radius of the corresponding nonrotating stellar model. 
This a is normalized to reduce to the standard definition 
used in perturbations of slowly rotating stars We 
must also define a generalization of the frequency of the r- 
mode. For a small-amplitude mode the time derivative of 



J22 is proportional to the frequency: J. 



(i) 
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we are led to define the following nonlinear generalization 
of the r-mode frequency 
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These expressions for a and uj are very stable numerically 
since they are expressed as integrals. jnj 

We have studied the growth of an r-mode as described 
above by solving numerically the nonlinear hydrody- 
namic equations for a rotating stellar model represented 
on a 64 x 128 x 128 cylindrical grid. We prepare initial 
data for this evolution by constructing first a rigidly ro- 
tating equilibrium stellar model. For simplicity we use 
the polytropic equation of state p = Kp 2 , which approx- 
imates the features of more realistic neutron-star mod- 
els reasonably well. The model most extensively stud- 
ied here is very rapidly rotating, with angular velocity 
= 0.635\/7rGpo where po is the average density of the 
nonrotating star of the same mass. The ratio of the ini- 
tial rotational kinetic energy to gravitational potential 
energy for this model is T"/|W| = 0.101. At the initial 
time we take the fluid density to be that of this equilib- 
rium stellar model, while the fluid velocity is taken to 
include a small amplitude r-mode perturbation: 



v = Qo0 - 



a R n ( — ) Be(Y&). 



(10) 



In our simulation we take the initial r-mode amplitude 
to be ao = 0.1. 

According to perturbation theory ||, the ratio of the 
GR growth time to the r-mode pulsation period is ex- 
pected to be tcr/Pt = 5.64 x 10 4 for the stellar model 
studied here. At this rate, growth of the r-mode is far 
too slow to be studied effectively via an explicit hydro- 
dynamic simulation. To overcome this, we artificially 
increase n m Eqs. (|) and (§) to a value such that 
tcr/Pt ~ 12.6, i. e., the strength of our GR reaction 
is about 4500 times stronger than it should be. This 
made the growth rate of the r-mode large enough so that 
it could be followed up to nonlinear amplitude with a 
reasonable amount of computing time, yet kept it slow 
relative to the dynamical and sound-crossing times. 
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FIG. 1. Evolution of the current quadrupole moment J22. 
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FIG. 2. Nonlinear evolution of the r-mode amplitude a. 



Figure |] shows the Re( J22) that results from the evolu- 
tion of this system. Time in these figures is given in units 
of the initial rotation period of the star: Pq = 27r/Oo. 
Figure [l] illustrates that the evolution of J22 is dominated 
by the sinusoidal r-mode oscillations. Figure || illustrates 
the evolution of the r-mode amplitude in this simulation. 
We see that the growth is exponential (as predicted by 
perturbation theory) until a ~ 2. Then some nonlin- 
ear process limits the growth; a peaks at a = 3.35 and 
then decreases rapidly. The evolution of the r-mode fre- 
quency lo defined by Eq. (||) is illustrated in Fig. [| The 
evolution of u> is quite smooth when the amplitude of the 
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r-mode is large: a > 0.5. At early (t < 10Po) and at late 
(i > 28Po) times when the r-mode amplitude is small, we 
see that other modes also make noticeable contributions 
to J22, and hence to w. 
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FIG. 3. Evolution of the r-mode frequency. 

Several authors jllj have suggested that the nonlinear 
evolution of an unstable r-mode could generate differ- 
ential rotation. This differential rotation would amplify 
existing magnetic fields in the star, and these in turn 
might significantly affect the evolution of the r-mode. 
We have explored this possibility by monitoring the av- 
erage differential rotation Af2. We find that Af2 grows to 
Afi w 0.41f2 at time t » 28Po and then decreases. Here 
Cl is the ratio of the angular momentum to the moment 
of inertia of the star. But the average value of Af2/fi 
may be misleading. Figure [| illustrates the spatial de- 
pendence of the azimuthally averaged angular velocity 
il(w,z) — J lW<p/27r at the time t — 25.6P). We see 
that the differential rotation is confined mostly to a thin 
shell of material near the surface of the star, and is par- 
ticularly concentrated near each polar cap. The bulk of 
the material in the star remains fairly rigidly rotating. 
Thus it appears that the magnetic fields generated by 
this differential rotation may have strong effects locally 
(e.g., leading to the ejection of material at the surface or 
along the rotation axis), but may not affect the global 
behavior of the r-mode as much as if the differential ro- 
tation were distributed more uniformly. 




FIG. 4. Spatial dependence of the azimuthally averaged 
angular velocity Q(zu, z)/Q. Rotation axis is the left edge of 
the figure; equatorial plane is the bottom edge. 

Energy and angular momentum are removed from the 



star by GR emission from the r-mode according to the 
expression 0,0: 
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Figure [| illustrates the evolution of the total angular mo- 
mentum J, the total mass M, and the kinetic energy of 
the fluid T. We see that M is essentially unchanged, 
while J decreases by about 40%. The numerical evolu- 
tion of J agrees with Eq. ( 11 ) to within a few percent jl3) . 
Figure || reveals that the kinetic energy of the star, T, 
continues to decrease for about two rotation periods after 
the rate of emission of J into GR falls to zero. 




FIG. 5. Evolution of the total angular momentum J/ Jo 
(solid curve), mass M/Mo (dashed), and kinetic energy T/Tq 
(dot-dashed) of the star. 

What nonlinear process is responsible for limiting the 
growth of the r-mode? Figure || reveals that the kinetic 
energy of the star, T, continues to decrease even after the 
rate of emission of J into GR falls to zero. This implies 
that the energy stored in the r-mode is not being radiated 
away as GR or being transferred to other macroscopic 
modes. For if it were, T would be more or less conserved 
once the GR losses become small. Rather we see that 
at t m 26Po, as the amplitude of the r-mode peaks, the 
predicted evolution of the total energy E due to GR loss 
(the dashed curve in Fig. ^|) significantly diverges from 
the numerical evolution of E. This nonconservation of 
E is due to the formation of shocks associated with the 
breaking of surface waves as illustrated in Fig. fj] 




FIG. 6. Evolution of the total energy E divided by the 
initial rotational kinetic energy To. Solid curve is the actual 
numerical evolution; dashed curve is the predicted evolution 
due to GR losses alone using Eq. ([H]) . 
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In summary then, we find that an r-mode can grow to 
relatively large amplitude under the influence of the GR 
reaction force. The hydrodynamic mechanism that acts 
to suppress the r-mode in our simulation is the forma- 
tion of shocks near the surface of the star. Since the GR 
reaction force in our simulation is too strong by a factor 
of 4500, it is still possible that slower hydrodynamic pro- 
cesses (like the transfer of energy to other modes) could 
limit the r-mode at smaller values of a. It is also pos- 
sible that the coupling of the r-mode to g-modes in real 
neutron-star matter (but absent from our barotropic sim- 
ulation) could also limit the growth at smaller a. How- 
ever if shock formation turns out to be the dominant 
suppression mechanism, then we expect the peak ampli- 
tude to be relatively insensitive to the strength of the GR 
coupling, iflojl In this case the dimensionless amplitude of 
the mode will peak at about a — 3.4. This implies that 
the spindown the star by GR will occur in about one 
tenth (oc 1/ o^ax) the time previously estimated || . 




FIG. 7. Density contours (at 10 



-n/2 



with n = 1, 2, 



in selected meridional planes at times t = 28Po and 29Po 
illustrate the breaking of surface waves. Shocks at the leading 
edges of these waves appear to be the primary mechanism that 
suppresses the r-mode. 

During our simulation about 40% of the initial angu- 
lar momentum and 50% of the initial rotational kinetic 
energy is radiated away as GR. Most of this energy is 
radiated in a much narrower frequency band A/ w 0.05/ 
than had been expected Q. Further, the frequency of 
the r-mode is about 20% smaller than that predicted by 
simple perturbation theory. Both of these effects tend 
to make the GR emitted by this process more easily de- 
tectable by LIGO. Even at the end of our simulation the 
star is still rather rapidly rotating, and we presume that 
GR will again drive the r-mode unstable and a second 
episode of spindown will occur. During our simulation 
about 16% of the initial rotational kinetic energy of the 
star is dissipated by the shocks. This energy will be con- 
verted to heat in a real neutron star and this would delay 
cooling and the formation of a crust. A detailed thermal 
analysis will have to be carried out to determine exactly 
what other effects this thermal energy may have. 
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